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ABSTRACT 



The abundance of both amorphous and crystalhne sihcates in very small grains is 
limited by the fact that the Wfim silicate emission feature is not detected in the diffuse 
ISM. On the basis of the observed IR emission spectrum for the diffuse ISM, the observed 
ultraviolet extinction curve, and the 10;um silicate absorption profile, we obtain upper 
limits on the abundances of ultrasmall (a < 15 A) amorphous and crystalline silicate 
■ grains. Contrary to previous work, as much as ~20% of interstellar Si could be in 

a < IsA silicate grains without violating observational constraints. Not more than 
~5% of the Si can be in crystalline silicates (of any size). 
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a • 1. Introduction 



The idea of interstellar silicate grains can be traced back to the 1960's when Kamijo (1963) first 



' proposed that Si02, condensed in the atmospheres of cool stars and blown out into the interstellar 

^ space, could provide condensation cores for the formation of "dirty ices". It was later shown by 

Oilman (1969) that grains around oxygen-rich cool giants are mainly silicates such as Al2Si03 and 
Mg2Si04. Silicates were first detected in emission in M stars (Woolf & Ney 1969; Knacke et al. 
1969), in the Trapezium region of the Orion Nebula (Stein & Gillett 1969), and in comet Bennett 
19691 (Maas, Ney, & Woolf 1970); in absorption toward the Galactic Center (Hackwell, Oehrz, & 
Woolf 1970), and toward the Becklin-Neugebauer object and Kleinmann-Low Nebula (Gillett &: 
Forrest 1973). Silicates are now known to be ubiquitous, seen in interstellar clouds, circumstellar 
disks around young stellar objects (YSOs), main-sequence stars and evolved stars, in HII regions, 
and in interplanetary and cometary dust. 

The observed broad, featureless profile of the 10//m Si-0 stretching mode indicated that in- 
terstellar silicates are predominantly amorphous. With the advent of Infrared Space Observatory 
(ISO) and the improved ground-based spectrometers, the presence of crystalline silicates has been 
revealed in six comets (see Manner 1999 for a summary), in dust disks around main-sequence stars 
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(see Artymowicz 2000 for a summary), YSOs, and evolved stars (see Waelkens, Malfait, & Waters 
2000 for a summary), in interplanetary dust particles (IDPs) (Bradley et al. 1999), and probably 
also in the Orion Nebula (Cesarsky et al. 2000). 

If silicate grains with radii a < 15 A - which we will refer to as "ultrasmall" arc present, 
single-photon heating by starlight (Greenberg 1968) will cause these grains to radiate in the lOfim 
feature (see, e.g., Draine & Anderson 1985). The absence of a Wfim emission feature in IRAS 
spectra led Desert et al. (1986) to conclude that not more than 1% of Si could be in a < 15 A 
silicate grains, and on this basis recent grain models (Duley, Jones, &; Williams 1989; Desert, 
Boulanger, k Puget 1990; Siebenmorgen k Kriigel 1992; Mathis 1996; Li & Greenberg 1997; Dwek 
et al. 1997; Weingartner & Draine 2001; Li & Draine 2001) have excluded ultrasmall silicate grains 
as a significant grain component. Nondetection of the 10/um silicate emission feature in the diffuse 
ISM by ISO (Mattila et al. 1996) and the Infrared Telescope in Space (IRTS) (Onaka et al. 1996) 
appeared to confirm this conclusion. 

However, the presence of ultrasmall silicate grains can not be ruled out since the 10/xm silicate 
emission feature may be hidden by the dominant PAH features. In this Letter we place quantitative 
upper limits on the abundances of both amorphous and crystalline ultrasmall silicate grains by 
calculating spectra for such tiny grains heated by starlight, and comparing to measurements of 

the infrared (IR) emission of the diffuse ISM by IRTS (Onaka et al. 1996) and by the Diffuse 
Infrared Background Experiment (DIRBE) instrument on the Cosmic Background Explorer (COBE) 
satellite. The interstellar extinction curve and the 10//m silicate absorption profile are also invoked 
to provide further constraints. 

In §2 we discuss the optical and thermal properties of crystalline silicate material. In §3 we 
calculate the IR emission spectra of ultrasmall silicate grains and deduce upper limits on their 
abundances based on the constraints provided by the observed IR emission spectrum (§3.1), by 
the interstellar extinction curve (§3.2), and by the 10/xm silicate absorption profile (§3.3). In §4 
we discuss the effects of dust size distributions and the sources of uncertainty. §5 summarizes our 
main conclusions. 

2. Optical and Thermal Properties of Amorphous and Crystalline Silicates 

We approximate the grains as spherical and use Mie theory to calculate their absorption and 
scattering properties. The complex refractive index m(A) = m'(A) + i m"(A) of amorphous silicates 
is taken from Draine & Lee (1984). 

For crystalline silicates, we proceed as follows: for A < 0.3/um, we take m" from Huffman & 
Stapp (1973) for crystalline olivine (Mg,Fe)2Si04; for 0.3 < A < 6/Ltm, we take m" from "astronom- 
ical silicate" (Draine & Lee 1984); for 7 < A < 200/xm, we take m" from Mukai & Koike (1990) 
for crystalline olivine Mgi.8Feo.2Si04; for A > 200/xm, we assume m"(A) = m"(200/xm)(200/xm/A)^. 
The real part m'(A) is obtained from m"(A) using the Kramers-Kronig relation. 
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The specific heat of amorphoTis silicates is taken from Draine &: Li (2001). For crystalline 
silicates, following Draine &: Li (2001), we adopt a three-dimensional Debye model with B ~ 720K 
which closely reproduces the specific heat of crystalline Mg2Si04 measured by Robie, Hemingway, 
k Takei (1982). 

3. Upper Limits on Ultrasmall Silicate Grains 

3.1. IR Emission Spectra 

Wc will now derive upper limits on the abundances of ultrasmall silicate grains in the dif- 
fuse ISM based on comparison of the observed IR emission spectrum with the calculated emission 
spectrum of ultrasmall silicate grains. For this purpose, we require a line-of-sight with: 1) a good 
quality mid-IR emission spectrum; and estimates of 2) the starlight intensity xmmp (relative to the 
Mathis, Mezger, Sz Panagia [1983] radiation field), and 3) the hydrogen column density A^h- 

The diffuse ISM toward I = 44°20', b = -0°20' has been observed by DIRBE (Hauser et 
al. 1998) and the Mid-Infrared Spectrograph (MIRS) on IRTS (Onaka et al. 1996) has obtained 
the 4.7-11.7/Ltm spectrum. Li Sz Draine (2001) have modelled the IR emission from the dust in 
this direction. In Figure 1 we show the observed photometry and mid-IR spectrum together with 
the Li & Draine (2001) model, with a total gas column A^h = 4.3 x 10^^ cm~^, and a mixture of 
carbonaceous grains and silicate grains heated by starlight with XMMP ~ 2. This model is seen 
to be in good agreement with the PAH emission features seen by IRTS, as well as the DIRBE 
photometry. The total power radiated constrains the product xmmp x -^h, while the ~100/im 
wavelength of peak emission constrains Xmmp- 

Lacking a priori expectations concerning the form of the size distribution for ultrasmall silicate 
grains, we first consider grains of a single size a (radius of an equal- volume sphere). Suppose a 
fraction Ysi(a) of the total silicon abundance (Si/H 36 x 10~^; Grevesse & Sauval 1998) is 
contained in amorphous silicate grains of radius a. Applying the thermal-discrete method developed 
in Draine & Li (2001), we calculate the energy probability distribution function and the resulting 
IR emission spectrum. The resulting emission spectra are shown in Figure 2, where for each size a 
the value of Ygj has been chosen such that this abundance of ultrasmall amorphous silicate grains 
would not - by itself - exceed the observed emission at any wavelength. Grains with a < 10 A are 
heated by energetic photons to "temperatures" high enough to emit in the 10/xm silicate feature, 
and their abundances are thus limited by the IRTS spectrum. For grains with a > 10 A the emission 
peak shifts to longer wavelength, and the DIRBE 25/xm data provides stronger constraints. 

From Figure 2a it can be seen that, for a < 10 A, even if a fraction Ygi = 15% were in the form 
of amorphous silicates, the resulting 10/im Si-0 emission feature would not have been apparent in 
the IRTS spectrum. For a = lOA, Ygi < 25% is permitted. For 10 < a < 20 A, the DIRBE 25/Lim 
photometry allows up to Y^i < 30%, and for 20 < a < 25 A, Y^i < 40%. 
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Fig. 1. — Spectrum toward I = 44°20', b = -0°20' (see text). Diamonds: DIRBE photometry Heavy solid curve: 
MIRS/IRTS spectrum. Light solid curve: emission calculated for the carbonaceous/silicate grain model of Li & 
Draine (2001) heated by diffuse starlight with an intensity xmmp ~ 2 relative to the Mathis, Mezger, & Panagia 
(1983) radiation field. Reddened starlight is labelled "stellar"; curves labelled Bcarb, Bsii are contributions from "big" 
(a > 250 A) carbonaceous and silicate grains; Scarb, Ssii are contributions from "small" (a < 250 A) carbonaceous 
and silicate grains. 

Similarly, we can estimate Ysi(a) for ultrasmall crystalline silicate grains. As shown in Figure 
2b, crystalline Mgi.8Feo.2Si04 has pronounced features at 10.0, 11.1, 16.3, 18.6, 23.0 and 33.6/um. 
For a < 10 A, nondetection of the 10.0/um feature implies Ysi < 7%. For 10 < a < 20 A, an 
upper hmit Y^i < 15% is imposed by the DIRBE 25//m data, and for 20 < a < 25 A, Ysi < 20%. 
Interstellar crystalline silicates could presumably differ from the crystalline olivine Mgi.8Feo.2Si04 
considered here, resulting in changes in the strengths and positions of the infrared resonances (e.g., 
see Jager et al. 1998 for different crystalline silicate samples). Observations of comet Hale-Bopp 
(see Manner 1999) showed an emission feature at 11.2 rather than ll.l^m, and an emission feature 
at 10.0/im, evident in some spectra, was not always apparent. Our limits on a < 10 A crystalline 
silicates relied on the strength and position of the 10.0/im feature; since this feature in interstellar 
crystalline silicates might be either weaker or shifted in wavelength, it is appropriate to weaken 
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our upper limits for a < 10 A to Ygj < 10%. Below we will find that the lO^m absorption profile 
provides a more restrictive limit on crystalline silicates. 




X (yu,m) 

Fig. 2. — Contribution to the IR emission toward (44° < / < 44° 40', -0°40' < 6 < 0°) by: (a) amorphous sihcate 
grains with radii a=3.5A (Y^i = 15%), SA (Y^i = 15%), 7.5A (Y^i = 15%), lOA (Y^i = 25%), 12.5A (Y^i = 30%), ISA 
(Ysi = 30%), 20A (Ysi = 30%), and 25A (Y^i = 40%); (b) crystalUne silicate grains with radii a=3.5A (Yai = 7%), 
5A (Y,i = 7%), 7.5A (Y,i = 7%), lOA (Y^i = 15%), 12.5A (Y,; = 15%), 15A (Y,; = 15%), 20A (Y^i = 15%), 25A 
(Ysi — 20%). Diamonds: DIRBE photometry. Heavy solid line: 5-12/im spectrum observed by IRTS. 



3.2. Ultraviolet Extinction Curve 

Ultrasmall silicate grains make a steeply-rising contribution to the ultraviolet (UV) extinction. 
Weingartner & Draine (2001) and Li & Draine (2001) have obtained a carbonaceous-silicate grain 
model which reproduces both the observed extinction and the observed infrared emission. This 
grain model has only a small abundance Ygi ~ 1.3% of ultrasmall (a < 15 A) silicate grains. 
In Figure 3 we show the extinction obtained by adding ultrasmall silicate grains to this model. 
Since the ultrasmall grains are in the Rayleigh limit, the added extinction depends only on AYgi 
(independent of grain size), and applies to both crystalline and amorphous silicates since their 
optical properties at A^^ > 5//m^^ are taken to be similar (see §2). 
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As expected, the contribution to the A~ > 7fim region would be considerable for AYgi > 
20% (although the effects on the optical and near-UV extinction are negligible). In view of (1) 
the uncertainties in the observed interstellar extinction in this wavelength region (for example, the 
determinations of Mathis [1990] and Fitzpatrick [1999] differ from one other by ~ 5%), (2) the fact 
that the size distributions of other components (carbonaceous grains and a > ISA silicates) could 
be adjusted to reduce their contribution to the far-UV rise, and (3) the possibility that a mixture of 
ultrasmall silicates (both neutral and charged) might have a more gradual onset of UV absorption 
than our adopted silicate dielectric function (which has a relatively sharp increase in absorption 
near 6.5-7.0/im~^), it seems clear that the observed UV extinction does not rule out AYsi as large 
as ~ 20% on lines of sight with Ry = 3.1. 




A-i (Mm-i) 

Fig. 3. — Comparison of the average Galactic interstellar extinction curve (open circles: Mathis 1990; dotted 
line: Fitzpatrick 1999) with theoretical extinction curves (solid lines) obtained by adding the contributions from 
ultrasmall silicate grains with AYsi=5%, 10%, 15%, 20% and 25% of the solar silicon abundance (thin solid lines) to 
our best-fitting model (with Yai = 1.3%; see Weingartner & Draine 2001; Li & Draine 2001). Also plotted are the 
contributions (see Li & Draine 2001) from "Bsh" (a > 250A silicate); "Ssii" (a < 250A silicate); "Bcarb" (a > 250A 
carbonaceous); "Scarb" (« < 250A carbonaceous, including PAHs). 
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3.3. Silicate Absorption Profiles 

The lO/Liin silicate absorption feature provides another constraint on the abundance of crys- 
talhne sihcate grains of any size. A recent high signal-to-noise spectroscopic study of lines of sight 
towards Cyg 0B2 # 12 (diffuse ISM), Elias 16 (quiescent molecular cloud) and Elias 18 (embedded 
YSO) in the Taurus dark cloud appears to rule out any narrow absorption features between 8.4 
and 12.5/im (Bowey et al. 1998). Figure 4 plots the lOfim silicate absorption profiles resulting 
from addition of various amounts (AYsi) of ultrasmall crystalline silicates to our best-fitting model 
(see Weingartner & Draine 2001; Li & Draine 2001). The absence of narrow features near 10.0 
and 11.1/xm in the observed 8.4-12. 5jLtm extinction (Bowey et al. 1998) suggests an upper limit 
AYsi < 3% for crystalline silicates. Note that this result is insensitive to size as the grains are in 
the Rayleigh limit. 

As noted above, the precise wavelengths and strengths of IR resonance features in crystalline 
silicates vary with composition, and there could be more than one crystalline form present in the 
interstellar medium. Nevertheless, the lack of fine structure in the observed broad 10/im absorption 
feature makes it unlikely that more than 5% of the interstellar Si could be in crystalline silicates. 
This limit on crystalline silicates is more stringent than the limit on ultrasmall crystalline silicates 
obtained in §3.1. 

4. Discussion 

The upper limits on the abundances of ultrasmall silicate grains estimated from the IR emission 
spectra in §3.1 are for grains of a single size. Adopting a log-normal size distribution dn/dlna oc 
exp{— [ln(a/ao)]^/(2(T^)}, we have calculated emission spectra for various qq cind a. For this dis- 
tribution, a^dn/d\na peaks at ap = aoe'^'^^. For each distribution, we determined the maximum 
value of Ysi such that the infrared emission did not exceed either the IRTS 4.7-11.7;um spectrum 
or the DIRBE photometry. Results are shown in Table 1. It can be seen that Ygi increases with 
Up-, for 3 < ao < 6 A and 0.3 < ct < 0.6, Y^i varies in the range of 20%-55%. We take Y^i < 20% 
to be the limit on a < 15 A silicate grains. 

One uncertainty in our upper limit on ultrasmall amorphous silicates arises from uncertainty in 
the absorption cross section for 1 < < 4/xm~^ where the interstellar starlight spectrum peaks. 
We have assumed the ultrasmall silicate grains to be characterized by the refractive index estimated 
for amorphous "astronomical silicates" (Draine &: Lee 1984; Weingartner & Draine 2001). We note 
that the visual absorptivity of crystalline olivine material is sensitive to its Fe content. If - as 
suggested by the low Fe content observed in some circumstellar environments (see Waelkens et al. 
2000) and in comet Hale-Bopp (Crovisier et al. 1997) - crystalline astronomical silicate material is 
less absorbing than "astronomical silicate" (Draine & Lee 1984), then photon absorptions will be 
less frequent and the IR emission correspondingly reduced, in which case the infrared emission limit 
on ultrasmall silicates would be relaxed. However, for crystalline silicates, the strongest constraint 
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Fig. 4. — The lO/im silicate absorption profile (thick solid line) for the carbonaceous-amorphous silicate model of 
Li & Draine (2001) as well as profiles with an additional contribution from crystalline silicate grains with AYsi=2%, 
4%, 6%, 8% and 10% of the solar silicon abundance (thin solid lines). The narrow feature at ll.B^im is due to a 
lattice mode of crystalline graphite (Draine 1984). Dot-dashed line is the absorption due to amorphous silicate (i.e., 
with the contributions of other dust components subtracted). Also shown is absorption which would be contributed 
by crystalline silicate grains for various Y^i- Also plotted is the observed extinction profile for dust toward Cyg OB2 
#12 (Bowey et al. 1998) divided by iVn = 2 x 10^^ cm'^ (Whittet et al. 1997). 

comes from lack of sharp absorption features near 10 and llfim, so our upper limit is not affected 
by uncertainty in the l~4/im~^ absorption. 

The observed 3.3, 6.2, 7.7, 8.6, and 11.3//m emission bands (Mattila et al. 1996; Onaka et al. 
1996) require a significant fraction of carbon in polycyclic aromatic hydrocarbon molecules with 
a < I5A (the model of Li & Draine 2001 requires ~ 40 carbon atoms per 10^ H, or ^^10% of 
the total C abundance). In contrast. Desert et al. (1986) argued that not more than -^1% of the 
Si could be in silicate particles with radii a < 15 A. Such a paucity of ultrasmall silicate grains 
would contrast sharply with the apparent large abundance of small carbonaceous particles. The 
present study, however, finds that as much as 20% of the Si could be locked up in ultrasmall silicate 
particles without violating any observational constraints. 

We suspect that the analysis by Desert et al. (based on 2-13;um observations of M82 and 
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Table 1: Upper limits on the abundance Ygi of ultrasmall amorphous silicate grains placed by the 
absence of the 10/xm silicate emission feature (see §3.1, §3.3). 
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0.3 


3.9A 


20% 


4A 


0.3 
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7.9A 


20% 
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0.4 
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9.7A 


30% 


3A 
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6.4A 


20% 


4A 


0.5 


8.5A 


25% 


5A 


0.5 


10. 6A 


35% 


6A 


0.5 


12. 7A 


45% 


3A 


0.6 


8.8A 


30% 


4A 


0.6 


11. 8A 


40% 


5A 


0.6 


14.7A 


50% 


6A 


0.6 


I7.7A 


55% 



NGC 2023) underestimated the fraction of the absorbed starlight energy which is reradiated in the 
6-12/im PAH features. The present study is based on a line of sight with photometry extending to 
the far-infrared (see Figure 1). 

5. Summary 

Upper limits on the abundances of ultrasmall amorphous and crystalline silicate grains are 
estimated based on observed IR emission spectra (§3.1), the interstellar extinction curve (§3.2), 
and the 10/xm silicate absorption profile (§3.3). Our principal results are: 

1. For amorphous silicate grains, non-detection of the lOfim silicate emission feature imposes an 
upper limit on the fraction Ysi of interstellar Si in very small silicate grains: Ysi < 15% for 
radii a < 10 A, and Y^i < 20% for a = lOA, Ygi < 30% for 10 < a < 20 A, and Y^i < 40% for 
20 < a < 25 A. For a distribution of sizes, Ygi < 20% in grains with a < 15A. These limits 
are much higher than the earlier estimate by Desert et al. (1986) of Ygi < 1% for a < 15 A. 

2. The observed UV extinction curve in diffuse regions may allow as much as Ygi < 20% in 
ultrasmall silicate grains. 

3. For crystalline silicate grains, the absence of fine structure in the observed lO/im silicate 
absorption profile provides the strongest constraint, limiting the abundances of a < 1/xm 
crystalline silicate grains to Ygi < 5%. 

Toward high galactic latitude "cirrus" both the radiation field and the hydrogen column density 
are known, but mid-IR spectroscopy is generally unavailable. In the near future, the Space Infrared 
Telescope Facility (SIRTF) may be able to obtain 5-40/im spectra of cirrus regions to cither detect 
silicate emission or obtain more rigorous upper limits on the abundances of very small silicate 
grains. 
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